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We present a detailed study of the vibrational properties of Single Wall Carbon Nanotubes 
(SWNTs). The phonon dispersions of SWNTs are strongly shaped by the effects of electron-phonon 
coupling. We analyze the separate contributions of curvature and confinement. Confinement plays 
a major role in modifying SWNT phonons and is often more relevant than curvature. Due to their 
one-dimensional character, metallic tubes are expected to undergo Peierls distortions (PD) at T=OK. 
At finite temperature, PD are no longer present, but phonons with atomic displacements similar to 
those of the PD are affected by strong Kohn anomalies (KA). We investigate by Density Functional 
Theory (DFT) KA and PD in metallic SWNTs with diameters up to 3 nm, in the electronic temper- 
ature range from 4K to 3000 K. We then derive a set of simple formulas accounting for all the DFT 
results. Finally, we prove that the static approach, commonly used for the evaluation of phonon 
frequencies in solids, fails because of the SWNTs reduced dimensionality. The correct description of 
KA in metallic SWNTs can be obtained only by using a dynamical approach, beyond the adiabatic 
Born-Oppenheimer approximation, by taking into account non-adiabatic contributions. Dynamic 
effects induce significant changes in the occurrence and shape of Kohn anomalies. We show that the 
SWNT Raman G peak can only be interpreted considering the combined dynamic, curvature and 
confinement effects. We assign the and G~ peaks of metallic SWNTs to TO (circumferential) 
and LO (axial) modes, respectively, the opposite of semiconducting SWNTs. 
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I. INTRODUCTION 

Carbon nanotubes are at the center of nanotcclinol- 
ogy researcbii^. The determination of their structure, 
phonon dispersions, and Raman spectra is a most intense 
area of investigation since their discoverji^iii^i^. Single 
wall nanotubes (SWNTs) can be described either as gi- 
ant molecules or one-dimensional crystals. Due to their 
reduced dimensionality, confinement effects play a fun- 
damental role in shaping their physical properties, such 
as the metallic or semiconducting character or the elec- 
tron transport mechanisms. Optical phonons of SWNTs 
are extremely important. Indeed, they contribute the 
most intense features in the Raman spectra. These are 
commonly used to sort metallic from semiconducting 
SWNTs. Electron scattering by optical phonons also sets 
the ultimate limit of high field ballistic transport^i^iiS, 
due to hot phonon generatio n^°i^^ . Electron-phonon cou- 
pling (EPC) is the key physical parameter necessary 
to quantify the phonons interaction with electrons. In 
metallic SWNTs the EPC strongly affects the phonon fre- 
quencies, giving rise to Kohn Anomalies ( K A ) ^ iii^ v]J A ^ ;lg, 
and Peierls Distortions ^P'D)15^Ld»^9j2^L22^m.. A cor- 
rect understanding and a quantitatively precise descrip- 
tion of electron-phonon coupling, Kohn anomalies, and 
Peierls distortions in SWNTs is then of prime scientific 
and technological interest. 

Experimentally, the most commonly studied SWNTs 
have diameters from 0.8 nm, up to 3 nm, with hun- 
dreds to thousands atoms per unit cell. Covering such 



diameter range pushes the limits of computational ap- 
proaches. SWNTs optical phonons were calculated us- 
ing phonon zone foldingi^i^i^^i^l, tight binding^^i^i^i^ 
(TB), density functional theory (DFT)^ii^a9i20^i32^ 
and symmetry-adapted modelsSii^'^i^'^. The simplest 
approach to get SWNT phonons, consists in folding those 
of graphene^i^ii^i^. This requires the prior calculation 
of graphene phonons either by TB^iS6,38 ^ DFTi^J^i21i39 
or by other approaches^. Phonon zone folding (PZF) 
is not limited by the size of the SWNT unit cell. How- 
ever, even if SWNT electrons are well described by folded 
graphene^, this is not always true for phonons. We have 
shown that graphene has two KA at T and Ki^. Due 
to their reduced dimensionality, metallic SWNTs have 
much stronger KA than graphen o^^'^'^i^^ . Thus, folded 
graphene cannot give a reliable phonon dispersion of 
metallic SWNTs close to the Kohn anomalies, no matter 
how accurate the graphene phonons^^''^^. On the other 
hand, semiconducting SWNTs cannot have KA. So, even 
in this case, folded graphene may not correctly describe 
the phonons. 

These problems can be overcome by means of TB or 
DFT calculations. Indeed, the EPC effects can be taken 
into account in TB models^'i^!^. However, the EPC 
strongly depend on the TB parameterisation, and con- 
trasting values are found in literature^iiSi^i^i^. Also, 
most TB works failed to identify KA^^^S^. DFT di- 
rectly includes EPC effects, and, being an ab-initio tech- 
nique, does not rely on external adjustable parameters. 
However, this better accuracy is obtained at the price of 
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a significant increase of the computational time demand, 
which makes the evaluation of SWNT properties a chal- 
lenge. As a result, most DFT works focus on extremely 
small tubesi^i^ and/or achiral tubesiiii^i^l. Finally, in 
all the approaches used to compute phonon frequencies, 
the time dependent nature of phonons is neglected. This 
approximation, which results in a static description of the 
atomic vibration, is usually legitimate in tri-dimensional 
crystals^^i^, but, as we will show, does not necessarily 
hold for one-dimensional systems. 

In this paper we show that SWNTs are a major ex- 
ception, and that a correct description of their optical 
phonons can only be achieved by carefully considering 
the dynamical nature of phonons. We first analyze the 
static case and single out the effects of curvature and 
confinement. Confinement plays a major role in shaping 
SWNT phonons and is often more relevant than curva- 
ture. We present an electronic zone folding method al- 
lowing the static DFT calculation of confinement effects 
on phonon-dispersions and EPC of SWNTs of any diame- 
ter and chirality, and for any electronic temperature (Tc). 
We investigate KA and PD in metallic SWNTs with di- 
ameters up to 3 nm and for Te=4 to 3000 K. We present a 
simple analytic model exactly accounting for all the static 
DFT results. Finally, we prove that dynamic effects in- 
duce significant changes in KA occurrence and shape. We 
show that the SWNT Raman G peak can only be inter- 
preted considering the combined dynamic, curvature and 
confinement effects. We assign the 0+ and G~ peaks of 
metallic SWNTs to TO (circumferential) and LO (axial) 
modes, the opposite of semiconducting SWNTs. A fam- 
ily dependence of the G peak position in semiconducting 
tubes is observed. 

This paper is organized as follows. In Section II we 
introduce the background concepts and definitions. Sec- 
tion III presents our EPC calculations. Section IV and V 
present respectively the study of the effect of pure con- 
finement within the static and the dynamic approach. 
Section VI highlights the influence of curvature on the 
SWNT phonons. Finally, Section VII assigns the G+ 
and G~ Raman peaks of SWNTs. 



II. BACKGROUND 

Before presenting our numerical and analytic results, it 
is necessary to introduce and discuss a set of background 
concepts and definitions, which will be used thorough the 
paper, in order to highlight the differences between our 
approach and others in literature. 



A. SWNTs unit cell and graphene zone folding 

SWNTs are usually identified by means of their chi- 
ral indices n and m**'^. The determination of the chi- 
ral indices is done by unrolling the unit cell of the tube 
and comparing it to an infinite graphene sheet. The un- 



rolled unit cell is a rectangular sheet of graphene that 
can be described by means of the vectors Ch and T. Ch 
is known as the chiral vector and has the same length as 
the circumference of the tube. T is known as the trans- 
lation vector and defines the translational symmetry of 
the tube. The crystal lattice of graphene can be defined 
by two vectors, ai and a.2, forming a 7r/3 angle. Ch and 
T can be projected on ai and a2'*'^: 

(2m + n)ai {2n + m)a2 

Ch = nSLi + mSL2 T = , (1) 

"R "R 

where n, m are defined as the chiral indices and dp = 
gcd(2n -t- m, 2m + n) (gcd= greatest common divisor). 

Wavevectors in nanotubes have to be commensurate to 
the tube circumference. They can be represented on the 
unrolled unit cell of the tube, and arc in the form: 

k = k_L + k|| =/iKi+j/K2, (2) 

where K2 is parallel to T, its modulus K2 = 27r/T 
correspond to the length of the tube's Brillouin zone 
(BZ), and v can assume any real value; Ki is perpen- 
dicular to T, and /iKi must be commensurate to the 
SWNT circumference. Thus Ki = 27r/Ch, M G N, and 
/i G [Nc/2—1, Nc—1], Nc being the number of unit cells of 
graphene contained in the unrolled unit cell of the tube. 
The number of nodes of a wavefunction along the tube 
circumference is given hy i = 2/i^i^i^i^. 

Fig. [1] plots the wavevectors commensurate to a (4, 4) 
tube on the reciprocal space of graphene, showing that 
they are located on lines oriented along K2, spaced by Ki 
and including the F poinf^i^i^i^. For the generic {n, m) 
tube, we refer to this ensemble of cutting lines as ZF(^n,m) ■ 
The procedure of using wavevectors of graphene to de- 
scribe wavevectors of a SWNT is known as Zone Fold- 
ing (ZF). Using the correspondence set by ZF between 
graphene and tube wavevectors, it is possible to use the 
electrons and phonons of graphene to obtain the band 
structure or the phonon dispersions of SWNTsii^i^i^. 

For clarity, in the rest of this paper we will distinguish 
the electron and phonon wavevectors by labeling them 
respectively with the symbols k and q. 

B. LO-TO phonons in graphene, graphite and 
nanotubes 

The analysis of the SWNT Raman G band shows that 
the phonons deriving from the T~ E2g mode of graphene, 
of A symmetry in SWNTs, provide the dominant contri- 
bution to the G+ and G^ Raman peaks^i^. Thus, in 
this paper we focus on these modes and their phonon 
branches. 

In literature, these phonons are often labeled using dif- 
ferent criteria^i^ii^i^i^i^: (i) their symmetry^i^i^, (ii) 
the direction of the atomic displacements with respect to 
the SWNT axis, radius or circumference^^, (iii) the longi- 
tudinal or transversal character of the graphene phonon 



3 




FIG. 1: (color on-line) Three equivalent representations of the 
zone-folding lines for a (4,4) tube, Z_F(4 4), in the reciprocal 
space of graphene. Left: plot of Z_F(4.4) in the hep BZ of 
graphene; the primitive vectors bi and b2 are plotted in blue, 
the lines of ZF(4,4) are plotted in red. Center: plot of ^7^(4 4) 
in the hexagonal representation of the graphene BZ. On the 
right: plot of Z_F(4_4) emphasizing the vectors Ki and K2. 
The length of K2 corresponds to the length of the (4, 4) tube 
BZ. 

branch from which they originate^^, (iiii) the longitudinal 
or transversal character relative to the SWNT— . 

In this paper the following definition is used: a mode 
is longitudinal if e||q and transverse if e_Lq; where e is 
the phonon polarisation (i.e. the direction of the atomic 
displacements) and q its wavevector, which also defines 
the direction of the phonon propagation. 

In SWNTs, phonons can propagate only along the axis. 
As a result, modes with atomic displacements parallel to 
the axis are longitudinal, and those with atomic displace- 
ments perpendicular to the axis (i.e. along the circum- 
ference) are transversal. 

In achiral SWNTs, symmetry forces the displacements 
of modes originating from graphene T — E2g to be either 
parallel or perpendicular to the axis^i^. Thus, in arm- 
chair and zigzag SWNTs these are exactly LO and TO. 

In chiral SWNTs, phonon polarizations are expected 
to be influenced by chiral angle and C-C bonds ori- 
entation^^. Thus, the modes originating from the 
graphene T — E2g should have only a prevalently-LO or a 
prevalently-TO character. Indeed, Ref. [Sll proposed the 
LO/TO classification to be meaningless in chiral SWNTs. 
However, in Section VI we will show that for any SWNTs 
with diameter > 0.8 nm, these modes are almost exactly 
LO and TO. We thus consider the LO/TO classification 
meaningful and convenient, also because it avoids differ- 
ent labeling in chiral and achiral SWNTs. 

C. Kohn anomalies and Peierls distortions 

In a metal, for certain phonons with a wavevector 
connecting two points of the Fermi surface it is possi- 
ble to have an abrupt change of the electronic screen- 
ing of the atomic vibrations. This results in a sud- 
den softening of the phonon-frequencies, which is called 
Kohn anomaly52*54,^ metallic SWNTs the Fermi sur- 
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FIG. 2: Geometrical condition for the KA onset in metallic 
SWNTs: for q=0 a point of the Fermi surface is connected to 
itself. This corresponds to a KA at F. q — K connects the 
two distinct Fermi points K and K', giving a second KA. ep 
is the Fermi energy. We label the band with positive slope 
L, and that with negative slope R. In electron transport, 
this notation indicates bands where electrons enter the tube 
respectively from the left and right electrode. 



face consists of two points only. Thus, KA can occur 
only for phonons with q = or q connecting the two 
Fermi points (Fig. [21). In graphene, a necessary condition 
for KA is a significant non-zero EPC (between electrons 
near the Fermi energy) for phonons at q =r or q = K— . 
Due to their reduced dimensionality, metallic SWNTs 
are expected to have stronger KA than graphite and 
graphene ^^"i^i^ii^ii^i^, for the corresponding phonons. 

One dimensional metallic systems are predicted to be 
unstable at T=0 K^. In fact, in ID metals it is always 
possible to find a lattice distortion that (i) opens a gap 
in the electronic structure, and (ii) such that the energy 
gained by gap opening compensates the strain induced 
by the deformation. This is called Peierls distortion^^. 
The possible presence of PD in SWNTs was argued as 
soon as their metallic behavior was predicted^i^. Their 
geometry was investigated by several authors with con- 
trasting results. It was proposed that they follow the 
atomic displacement patterns of LA phonons^ii^, of Ra- 
dial Breathing Modes (RBM)^^, of sohtwistons^i, or of 
optical phonons with q = P or g = 2fcpHii^a2ii^aai2£, 
where fcp is the wave-vector for which the electron gap 
is null. Note that the condition q = 2kF is equivalent to 
q = K in Fig. [2 

For T > OK, due to thermal excitation, states above 
the Fermi energy are populated. Thus, increasing the 
temperature reduces the energy gained by gap open- 
ing. PD are then stable only below a critical temper- 
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ature TpD • This temperature is also known as the metal- 
semiconductor transition temperature, and its determi- 
nation is crucial to understand the interplay between the 
possible onset of superconductivity and PD — 

For T > TpD, the energy gain achieved by gap opening 
is not sufficient to compensate the elastic strain, so the 
lattice will not undergo a permanent distortion. How- 
ever, phonon modes having a pattern of atomic displace- 
ments corresponding to the PD can show an important 
softening. This effect is proportional to the energy gained 
by gap opening, and inversely proportional to the tem- 
perature. This mechanism, described by some authors 
as a Peierls-like mechanisntiiii^, is, in fact, a particular 
case of Kohn Anomaly^. 



D. Temperature Effects 

Let us call Tc and Tj the electronic and the ionic tem- 
perature. Ti corresponds to the energy associated to the 
atoms vibrations around their equilibrium positions. Tc 
fixes the electronic states population close to the Fermi 
energy ep: 



1 



= (ck,„-CF)/fcBT<, _|_ I' 



(3) 



where /k,n is the Fermi-Dirac occupation of the n"* elec- 
tronic band at k, and /cb is the Boltzmann's constant. 

In metals, a finite electronic temperature always 
changes the electronic states occupation, while in mate- 
rials with an electronic gap Eg, changes of the electronic 
occupancies are relevant only for k^T^ > Eg. 

The Peierls-distortion temperature Tpo depends only 
on the electronic structure of the system, and not on the 
lattice thermal energy. Thus, Pcicrls distortions are fea- 
tures of the phonon dispersions which directly depend on 
To. On the other hand, phonons can have a temperature 
dependence due to anharmonicity^^, which is determined 
by Ti. 

It is thus essential to distinguish the different effects 
of Te and Tj on the SWNT phonons. At thermal equi- 
librium, Ti = Tc — Ts, where Ts is the sample tempera- 
ture. Thus, the correct phonon frequencies can be recov- 
ered only by taking simultaneously into account the con- 
tributions of Kohn anomalies and anharmonicity. The 
calculation of the anharmonic effects will be discussed 
elsewhere^. To a first approximation, the anharmonic 
contribution is similar for phonons of similar energy. On 
the contrary, a key result of this paper will be to show 
that To selectively acts only on the TO and LO phonon 
branches in metallic SWNTs, see Section IV and V. 

DFT, in its original formulation, is a ground state the- 
ory, thus Tc = OK. DFT can be extended to include a fi- 
nite real T^. This is done using the Mermin potential^ 
and populating the electronic states with the Fermi-Dirac 
distribution, Eq. [3] Because of the presence of KA and 
PD, the phonon frequencies of metallic SWNTs depend 



on Tc and DFT calculation should always be done using 
the Mermin potential. 



E. The dynamical matrix 

We call Oq^*^ the dynamical matrix projected on the 
phonons normal coordinates. Within time dependent 
perturbation theory, for a SWNT this is defined as (see, 
e.g., Eqs. 4.17a, 4.23 in Ref.^II): 



r\ 2T y—^ f \Dtk+q)n,km\'^[fk,m — fk+q,n] J, 



BZ 



<^k,m — £k+q,n + f^q + »7 



J An*{r)K{r,r')Anq{r') drdr' 
J n(r)AV^(r) dr, 



(4) 



where 27r/T is the length of the tube BZ, „ is a sum 
on all the possible electronic transitions, J^^ dk is an 
integral over the one dimensional tube BZ; efc^„ is the 
energy of the electronic Bloch eigenstate with wavevec- 
tor k and band index n; 7 is a small real constant; n(r) 
is the charge density; K[r,r') = 8'^EYi^c[n]/5n{r)5n{r'), 
where -Ehxc['t^] is the Hartree and exchange-correlation 
functional, and A^l/^ is the second derivative of the 
bare (purely ionic) potential; D(^k+q)n.km is the electron- 
phonon coupling matrix element: 



D 



(k-{-q)n.km 



(k + q,n\AVq[Ang]\k, 



(5) 




where AVq and Auq are the derivatives of Kohn-Sham 
potential and charge density with respect to displace- 
ment along the phonon normal coordinate; \k,n) is the 
electronic Bloch eigenstate of wavevector k and band n. 
The phonon frequencies ujq are derived from O^ aSi^: 



(6) 



where M is the atomic mass of carbon. 

It is important to note that Eq. [4] introduces a depen- 
dence of the dynamical matrix on the term huq+ij. This 
is a direct consequence of the perturbative approach used 
to derive Eq. |4l where the phonon is described as a time- 
dependent perturbation of the system®^. This accounts 
for the dynamic nature of vibrations. In the rest of the 
paper, we will refer to the inclusion of the term hujq + ij 
in the expression of the dynamical matrix as "the inclu- 
sion of, or "the description of dynamic effects (DE). 
This definition is introduced in opposition to the static 
approximation, which is described in Section II-E.2. 



1. Dynamic effects 

Dynamic effects induce a correction in the denomi- 
nators of terms proportional to the EPC in Eq. 21 In 
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materials where the electronic gap is larger than the 
phonon energy, the dynamic effects are negligible. In 
metals, dynamic effects are negligible if the phonon en- 
ergy is smaller than the electronic smearing for which the 
phonon calculations converge. In general, this is true for 
tri-dimcnsional metals^i^. Thus, in this case the con- 
tribution of the dynamic effects to the dynamical matrix 
is negligible. That is why dynamic effects are not in- 
cluded in the most common ab-initio codes for phonon 
calculations^!^. However, we will show that this is not 
necessarily true for systems with reduced dimensional- 
ity. For example, as explained in the previous Section, 
one dimensional metallic systems undergo Peierls distor- 
tions, the phonon frequencies depend on the electronic 
temperature, and the phonon energy can be of the same 
order, or even larger than the electronic temperature. In 
this case, it is conceptually wrong to neglect the dynamic 
effects a priori. 

Since metallic SWNTs can undergo Peierls distortions, 
the role of the dynamic effects on phonon frequency has 
to be carefully considered and investigated. This will be 
done in Section V. 



2. The Static Approximation 

The Born-Oppenheimer— , or adiabatic, or static ap- 
proximation is equivalent to dropping the phonon energy 
hujq in Eq. |4pA. Implementation of the static approxi- 
mation results in what we call static DFPT. To the best 
of our knowledge, all existing DFPT codes implement 
static equations. Static DFPT gives the same results as 
other static approaches, such as frozen phonons^^, where 
phonons are obtained evaluating the energy associated to 
a static lattice distortion. The expression of Og within 
the static approximation is obtained from Eq. [3] by set- 
ting to zero the phonon energy ujq and the term ij: 



e - 



fk.rn fk+q,: 
FBZ ^k,m ^ ^k+q,n 



{k~\-q)n,km \ 



dk 



j An;(r)ii'(r,r')Ang(r') dvdr' 
j n(r)AV(r) dr. 



(7) 



Unlike in the dynamic case, within the static approxi- 
mation Qq is real. Thus, the phonon frequencies tOq are 
simply given by: 



(8) 



3. Non analytic terms 

For states close to the Fermi energy, the denominator 
in Eq. [7]goes to zero, resulting in the possible presence of 
singularities in the phonon dispersion, as we have shown 



in graphene and graphite^^. Similar singularities occur 
in the real part of Eq. 2] when the energy of the phonon 
matches an electronic transition (see, e.g., Eq. D.17 in 
Ref. . To investigate such singularities in SWNTs it is 
convenient to split the dynamical matrix into an analytic 
and a non-analytic component: 



^9 



(9) 



The non analytic terms in the dynamical matrix of 
metallic tubes can be obtained by restricting the integral 
in the first term of Eqs.[l]and[7]to a set of k-points very 
close to the Fermi vector, and by limiting the sum over 
the electronic bands, taking into account only transitions 
between the bands crossing the Fermi point. Labeling 
with L the electronic band with a positive slope in Fig. [21 
and R that with negative slope, from Eqs.|l]and[7]we get: 



2Ar,KT 



2tt 



E 



m^n—L^R 



\D{K+k' +q)n,(K+k')mY 



fK+k', 



f K+k' +q,n 



<^K+k' 



SK+k'+q,: 



huj„ 



-dk' 
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(10) 



if dynamic effects are included, k has a small but finite 
value. In this case, divergencies occur when eK+k',m — 

£K+k'+q,n ± fii^q — 0. 

Within the static approximation we get: 



2Ar/KT 



2tt 



E 



m.n—L.R 



fK+k'jn — .fK+k'+q,n 
-k ^K+k' ,m ~ ^K+k'+q,n 



\D(K+k+q)n,(K+k')m\^dk' . 



(11) 



Here divergencies occur when €K+k',m = £K+k'+q,n- 

For both Eq.fTOland Eq.[TTl divergences occur for g ~ 
OT q ^ K (Fig. (2). In the first case, Ar — 2, in the 
second case Ak = 1. Note that the distinction between 
the analytical and non-analytical part of is operative, 
not physical. In fact, Qq contains all the divergent terms 
of Qq but not only them. The cut between Qq and 0™ 
depends on k. 



4. Curvature and Confinement 

The differences between graphene and a SWNT can be 
described in terms of curvature and confinement. Cur- 
vature effects are due to the distortion of the C-C bonds 
in the SWNT geometry and in a change of the bond 
character with respect to planar graphene. Confinement 
effects stem from the reduced dimensionality and are due 
to the quantization of the electronic wavefunctions along 
the SWNT circumference. 

In this paper we concentrate on the description of the 
effects of confinement and curvature on the dynamical 
matrix and on the phonon frequencies of SWNTs. Con- 
sidering the effects of curvature and confinement as per- 
turbations on the dynamical matrix of graphene Qq, at 
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first-order: 

Gt = Og + Ocurv + Oconf , (12) 

where Qt is the dynamical matrix of a SWNT and Gcurv 
and Bconf account for curvature and confinement. The 
effects of pure confinement will be calculated in Section 
IV and V. The effects of curvature will be evaluated in 
Section VI. 





LO 


TO 


K 


\Dlr\'' = \Drl\'' 

\Dll\' = \Drr\^ 


2(Di^>F 















TABLE I: EPC \Df in SWNTs for phonons corresponding 
to the graphene r-E29 ( LO or TO ) and K-A'i phonons. 
Electronic states are L and R bands near the Fermi energy 
(see Fig. [Hi. {D^)f = 45.60 eV A'^ and {D^)f = 92.05 
eVA-2, from Refs.[l3li2- 



III. ELECTRON-PHONON COUPLING 

In Ref. lid we have shown that the EPC of SWNT of 
arbitrary chirality is weakly affected by the tube curva- 
ture, and that it can be computed by folding the graphene 
EPC. Neglecting curvature, the EPC of a (n, m) SWNT 
is related to that of graphene by the ratio of the unit 
cells: 



|D2||ChxT| = |2?2 



]V3 



(13) 



where D is the EPC of a (n, m) SWNT ; D is the 
graphene EPC; Ch x T the area of the unrolled tube 
unit cell; oq is the lattice parameter of graphene, and 
ag\/3/2 is the graphene unit cell area. This formula is a 
simple consequence of electron and phonon normalization 
in the two different unit cells. Its derivation is presented 
in Appendix A. 

Since T and Ch are perpendicular, |Cii| — nd, being d 
the tube diameter, from Eq. [H the EPC of a SWNT is: 



m = 



^0 



2TTTd' 



(14) 



Using Eq. [T3]in Eq . [TUl and [TTl it is possible to obtain 
the expression for the non-analytic part of the dynamical 
matrix within zone folding. Eg. 1101 becomes: 



V 



E 



\D 



(K-|-k'-|-q)n,(K+k')m| 



m^n—L^R 
fn+k'/m — /K+k'+q,i 



-fc eK+k',m — eK+k'+q,n + ^^q + ^7 

while Eg. [TT] becomes: 

m.n—L.R 
''^ /K+k',m — /K+k'+q,ri 



dfc', (15) 



^(K+k'+q)n,(K+k')mP 



-dk'. 



k eK+k',r?j — eK+k'+q,n 



(16) 



Note that, since in zone-folding the electronic states 
of the tube are mapped onto those of graphene, all the 
wavevectors in Eqs. I15I16I are now vectorial quantities, 
and not scalar as in Eqs. IIOIIII The integration is per- 
formed along the line of ZF(^^a.m) passing through K. 



Eqs. I15|16l also show that the dependence of 8q on the 
translational vector T disappears. This implies that, for 
metallic tubes, any dependence on tube chirality is lost as 
well, and all the information on tube geometry is given 
by the ^ dependence only. 

As already stated, the relevant EPCs correspond to 
the graphene phonons at T and K, calculated between 
electronic states L and R near the Fermi energy (Fig [2]). 
In Ref. [13 it was shown that the only optical phonons of 
graphene with a non-negligible EPC are the E2g phonon 
at r and the A[ at K. Thus, throughout the paper we 
will consider only the SWNT phonons corresponding to 
the graphene r-E2g and to the K-A'j^. For simplicity, 
we will label them T ( TO or LO ) and K. For these 
phonons, in Ref. [l3 it is shown that the EPC can be 
considered independent from k and assumes the values 
of Tab. H] 

The EPC of graphene (and, thus, of SWNTs) are de- 
fined in literature in different ways: (i) the derivative of 
the hopping integral with respect to the C-C distance 
{rj)^JlJ2Jl, (ii) the EPC matrix element D of Eq.[5l (in) 
D times the phonon characteristic length, (.g)^^. The 
values of 77, D, and g are related by the following expres- 

,.;„„^10.12. 



= 9/2772 and \g\^\D\ 



(17) 



The relation between 77 and D is obtained by means of 
a first neighbors tight-binding model^. 

A variety of contrasting EPC values are reported in 
literature. Tab |TT] compares literature results for the 
graphene EPC with those we experimentally measured 
and calculated in Refs. IT^ItoI . This shows significant 
discrepancies between previous works and experimental 
EPCs, which, on the contrary, are in excellent agreement 
with our DFT calculations. 



IV. CONFINEMENT EFFECTS WITHIN 
STATIC DFPT 

In this Section we present a zone folding technique for 
the calculation of pure confinement effects on the phonon 
frequencies of SWNTs. All the results and calculations 
performed in this Section are obtained within the static 
approximation 
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77 [eVA ^] 


1 7-1 l2 r -x r2 A — 2l 

\Dlr\ [eV^A ^] 


\gLRr [oV'A ^] 


Mahan^ 


4.8-6.2 


103-172 


0.092-0.15 


Parl4 


6 


164 


0.146 


Jiang^ 


1.94 


17 


0.015 


Perebeinos^ 


5.29 


126 


0.112 


Koswatta^i 


6.00 


162 


0.144 


DFT^"- 


4.57 


94 


0.083 


Expcrimentali^ 


4.50 


91 


0.081 



TABLE II: Comparison between calculated and experimental 
EPCs for the F LO phonon, in its equivalent rj, |Dp and \gf 
expressions. All data in literature were calculated with tight 
binding. The relation between rj and D is obtained by means 
of a first neighbors tight-binding model. 



A. Electron and Phonon Zone Folding 

A simple and effective method exists to calculate the 
SWNTs electronic structure taking into account confine- 
ment only. Indeed, in Section II-A, it was shown that it 
is possible to fold the wavevectors of graphcnc into those 
of a {n, m) SWNT. This technique can be used to obtain 
the electronic band structure of SWNTs from that of 
graphene^i^i^i^. We refer to this technique as Electron 
Zone Folding (EZF). Using EZF it is possible to obtain 
the complete electronic structure of any SWNT doing the 
calculation just on the 2 atoms unit cell of graphcnc. EZF 
relies on the exact correspondence between the unrolled 
SWNT unit cell and graphene, therefore all the effects 
of curvature are completely neglected. Despite the sim- 
plicity of this technique, for SWNTs with d > 0.8 nm, 
the electronic structure calculated using EZF does not 
significantly differ from that calculated by DFT on the 
SWNT unit cell'*°. 

The same procedure used to fold electrons can be used 
to fold phononsi i^i^^i'^^ . We refer to this as Phonon 
Zone Folding (PZF). As already mentioned, KA affect 
graphene and metallic SWNTs differently, and do not af- 
fect semiconducting SWNTs at all. Thus, PZF is not 
suitable for the description of the phonon dispersion of 
SWNTs close to the anomalies. 



B. Phonon calculation within EZF 

Even if PZF cannot be used for the description of 
the Kohn anomalies, this result can still be achieved 
by means of a folding approach. Eq. [7] shows that the 
dynamical matrix of a SWNT depends on its electronic 
structure through: (i) the electronic charge density and 
(ii) an integral of the electronic states energies over the 
BZ. Numerically, both these quantities are obtained by 
performing discrete sums on a set of k-points over the 
tube BZ. These sums can be replaced by sums on points 
of in the two-dimensional graphene BZ. The 

evaluation of Qq within EZF completely neglects the ef- 
fects of curvature. Thus, using Eq. [T^l we find that 



SWNT dynamical matrices computed within EZF are in 
the form: 

e^-* ^Qg + econf (18) 

In this Section, we use DFPT to calculate the SWNT 
dynamical matrices within EZF. We refer to these calcu- 
lations as EZF-DFT. 

A further increase in the speed of calculations can 
be obtained from an appropriate choice of the BZ sam- 
pling. The graphcnc Fermi surface consists of the two 
inequivalent points K and K'. If the ZF(^n,m) lines cross 
these points, the tube is metallic, otherwise it is semi- 
conducting. The TT and tt* bands of graphene near the 
Fermi points have a conic shape, so the electronic bands 
of metallic SWNTs close to K and K' are linear, and 
ek' = if3k' . This implies that for metallic tubes the inte- 
gral in Eq. [7] contains terms that diverge like 1/fc. These 
singularities are responsible for the KA onset, and have 
to be properly described. The most efficient way to nu- 
merically integrate 1/k divergences in one dimension is to 
use a logarithmic sampling^^. Thus, to optimize the effi- 
ciency of our calculations, we sample the EZF lines with 
evenly-spaced grids of k-points in regions away from the 
Fermi surface, and with denser logarithmic grids when 
crossing K or K'. Details on the k-point sampling used 
for EZF-DFT calculations are reported in Tablllll 

1. Comparison to other computational techniques 

It is important to discuss how our EZF-DFT re- 
sults compare to other computational approaches in lit- 
erature. These can be divided in the following cat- 
egories: PZF,^^ tight binding symmetry 
adapted methods ^^^^iZi and BFTU^^^I^o^I^^. 

The advantages of EZF-DFT over PZF have already 
been discussed in the previous Sections. The advantage 
of EZF-DFT over conventional DFT is the possibility of 
using extremely dense k-point grids and to converge cal- 
culations for a smaller Tc. This is evident from a survey 
of previous DFT calculations. We start examining the 
DFT approach used in Refs. ,17,18..27,51. . These calcu- 
lations were performed on the real tube unit cell. Thus 
the use of dense k-point grids to converge the phonon 
frequency at realistic was computationally extremely 
expensive. Calculations of Refs. I27ll5ll correspond to BZ 
sampling of at most 3 non equivalent k-points in the unit 
cell. The comparison between the number of k-points 
used in Refs. 27,51, and those used in this and other 
works, seriously challenges the convergence of the calcu- 
lations in Refs. [27ll5ll . 

Refs. [TtIITsI used much denser k-point sampling (up to 
57) and smeared electrons with an Hermite-Gauss func- 
tion at a fictitious^ electronic temperature Tf = 1160K. 
In principle, these calculations should reproduce the 
Tc = results, showing PD in metallic tubes. However, 
these authors observed no PD, but just a softening of the 
phonon frequencies. This happens because their chosen 
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Ref. 
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for which curvature effects are not neghgible. This point 
wiU be discussed in more depth in Section VI. 

EZF-DFT has also to be compared with tight bind- 
ing and symmetry adapted calculations. Non-ab-initio 
techniques, such as tight binding, offer another way to 
overcome the extremely large CPU time required by 
DFT. Most TB models are unable to describe KA and 
Pq28^29^^ a detailed study of EPC and effects on 
SWNT phonon within TB was reported by Barnett et 
alJ^iS. This study investigated KA and PD in SWNTs, 
but the results are in contrast with those obtained by 
DFT^°. Another possibility for the efficient calculations 
of phonons in SWNTs is the development of symme- 
try based force constants methods^ and lattice dynamic 
models based on the screw symmetry of SWNTs^li^^iIi. 
So far, only TB implementations of symmetry adapted 
techniques were applied to the calculation of the vibra- 
tional properties of SWNTs, which thus may lack the 
precision of DFT. The implementation of screw symme- 
tries into a DFT scheme would be highly desirable, but 
would require writing a new dedicated code. On the con- 
trary, another advantage of EZF-DFT is that it can be 
performed with the already existing packages. Finally, it 
is very important to stress that all the techniques pre- 
sented in this Section rely on the static approximation, 
and thus neglect dynamic effects. 



TABLE III: k-points sampling used in our EZF-DFT calcu- 
lations and data from other authors. The number of k-points 
rife refers to the number of non equivalent points the tube BZ. 
These Uk points fold into UkNc points in the graphene BZ, 
Nc being the number of unit cells of graphene present in the 
tube unit cell. The column Te is the electronic temperature 
used in the calculations (in Kelvin). Values marked with * 
correspond to calculations where electrons are smeared using 
a Hermite-Gauss instead of the Fermi-Dirac function. The 
M/S column specifies the metallic or semiconducting nature 
of the tubes. The In column specifies if a logarithmic grid is 
used. 



value of Tf is too large to reach the limit of To 0, 
and describes the system at a finite, but undefined, value 
of Te. Calculations at lower Te were prevented by the 
huge CPU time requirements. DFT calculations at lower 
values of Tg have been reported only for extremely small 
SWNTs in Refs. 19,20. Here, they calculated the phonon 
dispersion of SWNTs with d ~ 4A at To as low as 137K, 
showing the presence of KA and PD. The comparison be- 
tween the k-points sampling and the electronic tempera- 
ture in literature and those used in this paper is reported 
in Tab. UTTl 

As wc will show in the next Sections, the calculations 
of the phonon dispersion of large SWNTs at To as low as 
4K can be be easily done by means of EZF-DFT. This is 
achieved at the price of neglecting the effects of curvature, 
which prevents EZF-DFT to be used for the study of the 
radial breathing modes, and results in inaccuracies in the 
determination of the phonon frequencies of those tubes 



C. Numerical Results 

In this Section, we use EZF-DFT to perform a system- 
atic study of the phonon dispersion of SWNTs, paying 
particular attention to the description of KA and PD. 
Calculations are performed within the static approxima- 
tion. 



1. Phonon Dispersions 

Here we use EZF-DFT to calculate phonons of SWNTs 
of different diameters and chiralities, focusing on diame- 
ter and To dependence of the KA induced phonon soft- 
ening. 

We start with the comparison of phonon zone folding 
(PZF) and EZF-DFT for metallic tubes. Fig. [3] compares 
the phonons of a metallic (11,11) tube calculated using 
PZF with those obtained by EZF-DFT. To improve read- 
ability, we only plot phonon branches corresponding to 
the lines of ZF(ii,n) crossing T and K. AU PZF and EZF- 
DFT calculations are performed using DFPT^^. The pa- 
rameters for PZF calculations are the same as those used 
for the calculation of graphite in Ref. 

PZF calculations use the phonons of graphene, thus: 
(i) at r no LO-TO splitting is observed, and (ii) the 
kinks and the phonon softening caused by KA at T and 
K are the same as those calculated in Ref. [l^. Since the 
phonon frequencies of graphene converge at an electronic 
temperature To SOOOif , no temperature dependence 
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(11,11) — 



KA@r 
KA@K 





0.0 



0.0 0.2 0.4 0.6 0.8 1.0 

Phonon Wavevector {2k/T units) 

FIG. 3; (color on-line) Phonons of the (11,11) metallic SWNT 
calculated using (a) Phonon Zone Folding (PZF), (b) Elec- 
tronic Zone Folding (EZF-DFT) at Tc = 4K, (c) 315K, and 
(d) 3150K. Branches affected by (Red) KA at T, (Blue) KA 
at K. 



of the phonon frequencies calculated by PZF can be ob- 
served for 0K< To <3000K. EZF-DFT calculations show 
a number of fundamental differences with respect to PZF. 
Fig. [3] (b,c,d) show that the phonon softening induced 
by KA at T and K strongly depends on To, with the 
phonon frequencies decreasing for decreasing T^.. Most 
importantly, we observe that at T only the LO mode is 
affected by KA. This causes the splitting of the LO-TO 
modes. As already stressed, this splitting cannot be ob- 
served using PZF. Due to the temperature dependence of 
the LO mode, the LO-TO splitting increases for decreas- 
ing To- For Tc ~ 3150K, the KA effects are completely 
removed, and the LO-TO splitting is zero. Finally, in Fig. 
m we compare the whole phonon dispersion of the (11,11) 
tube computed using (a) PZF, (b) EZF-DFT at = 4K 
and (c) To = 315K. Fig. H shows that PZF and EZF- 
DFT differ only for the phonon branches affected by KA, 
and are indistinguishable for all other phonon branches. 
Thus, in the phonon dispersion of metallic SWNTs only 
the modes affected by KA depend on Te, and PZF can be 
safely used for the calculation of all other phonon modes. 

It is also interesting to compare PZF and EZF-DFT for 



0.2 0.4 0.0 0.2 0.4 0.0 0.2 0.4 

Phonon Wavevector {2n/T units) 



FIG. 4: (color on-line) Comparison of the complete phonon 
dispersion of a (11,11) tube calculated using (a) PZF, (b) 
EZF-DFT at n = 315K, (c) EZF-DFT at To = 4K. The 
three calculations differ only for the branches affected by KA. 



semiconducting tubes. Fig. [5] compares PZF and EZF- 
DFT for a semiconducting (19,0) tube. This has the same 
diameter (1.5 nm) as the metallic (11,11) in Fig. [31 As 
for FiglJl in Figl^lwe only plot selected phonon branches, 
in particular those corresponding to the line of ZTjig o) 
crossing T and those corresponding to the two lines that 
lay most closely to K. KA cannot be present in semi- 
conductors, thus no KA should be observed in the (19,0) 
tube. Since ZT(i9 g) does not include K, and the empty 
and occupied electronic states are separated by a gap 
Eg ^ 0.5 eV, no temperature effects can be observed for 
To < Eg/kB - 6000K. 

Fig. [5] shows that the discrepancies between PZF and 
EZF-DFT are extremely reduced with respect to those 
of the metallic (11, 11) tube, the only noticeable differ- 
ence being the absence of LO-TO splitting in the PZF 
calculations. 



D. LO-TO splitting 

Within EZF-DFT, the doubly degenerate T - £^29 
phonon of graphene splits for both metallic and semi- 
conducting tubes. This splitting is a direct consequence 
of confinement in SWNTs. 

Defining ^lo^to) ^^'^ phonon frequency of the 
LO (TO) mode at q=0 for a semiconducting (metal- 
lic) SWNT, we also observe that w^o > ^^lq, while 



> uj: 



TO- 



This is explained as follows. From Tab. U 
the EPCs between L and R bands at T have a finite 
value for LO modes, and are null for TO modes. Since 
the EPC associated with TO phonons is null, from Eq.fTT] 



we obtain that cj^q^'' 
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TO' 



where uj', 
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FIG. 5: (color on-line) Phonons of the (19,0) tube calculated 
using (a) PZF on a q-points grid sampling Z_F(i9_o) i (b) EZF- 
DFT. Branches that in graphene and metallic tubes would 
be affected by the KA anomaly at T are in red, while those 
closest to the KA at K in graphene are in blue. 



phonon frequency of the TO mode at T for graphene. On 
the other hand, LO modes have a finite EPC. Here it is 
important to remind that what determines the physical 
properties related to EPC is not the EPC per se, but the 
product of the EPC times the electronic density of states 
(DOS) (see also Appendix A). Close to the Fermi energy, 
for semiconducting SWNTs, the DOS is null, while for 
metallic SWNTs it is a constant, and for graphene is zero 
at Fermi energy (ep) and increases linearly for e > ep and 
e < ep. It follows that DOS% < DOSf/P^ < DOS^, 
where the superscripts S, M, and graph refer respec- 
tively to semiconducting SWNTs, metallic SWNTs, and 
graphene, and the subscript ep indicates that the DOS 
is evaluated for energies e close to the Fermi energy ep. 
From Eq. [Til and using the same notation as the TO 



graph 

LO 
graph 



Since 



it follows that 



LO- 



modes, we then obtain uj^Jq < uj' 
in graphene we have 

'^LO > ^TO and < '^TO- 

With EZF-DFT, we also find that the atomic displace- 
ments of the modes deriving from the graphene T — E2g 
are parallel or perpendicular to qy within an angle less 
than 3.25°. The effect of confinement is then, for any chi- 
rality, to split the F — i?2g mode of graphene into almost 
exactly LO and TO modes. 



1. Diameter, electronic temperature, chirality and family 
dependence of phonon softening 

We now consider the dependence of the phonon disper- 
sions on SWNT chirality, diameter and electronic tem- 
perature. 



EZF KA@r 
EZF KA@K 
EZF 




0.0 0.2 0.4 0.6 0.8 1.0 

Phonon wavevector {InU units) 

FIG. 6: (color on-line) Phonons of the (11,11), (18,0), and 
(15,6) tubes calculated using EZF-DFT. All tubes have d ~ 
1.5 nm. Calculations are done at Te = 315K. Phonon 
branches affected by KA at T are in red, those with KA at K 
are in blue. 



Fig. [B] plots the phonon dispersions at Te=315 K for 
the (11,11) armchair tube, the (18,0) zig-zag tube, and 
the (15,6) chiral tube. All these tubes are metallic, and 
have a similar diameter of 1.5 nm. Note that in SWNTs, 
defining dR, = gcd(2n-|-TO, 2TO-|-n), kp = if (n— m)/3dR 
is not an integer (e.g. the (18,0) tube), and kp = 1/3, 
in 27r/r units, if (n — TO)/3(iR is an integer, (e.g. (11,11) 
and (15,6) tubes)^. This explains the different KA po- 
sitions in Fig. [HI We observe that, even if the shapes of 
the phonon dispersions are different, the frequency of the 
phonons corresponding to the graphene T — E2g-LO and 
K— A'j^ are the same in all cases, and are significantly 
softened with respect to graphene. Thus, the phonon 
softening due to KA does not depend on chirality. 

On the other hand. Figs. I7I4I clearly show that KA- 
induced phonon-softening strongly depends on diameter 
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FIG. 7: (color on-line) Phonons of the (6,6), (11,11), and 
(18,18) SWNTs calculated using EZF-DFT. Tube diameters 
are 0.8 nm, 1.5 nm, and 2.4 nm respectively. Tc is 315K. 
Phonon branches affected by KA at F are in red, while those 
with KA at K are in blue. 



and Te. Fig. [7] compares the phonons for three differ- 
ent armchair tubes at the same Tc. It indicates that the 
smaller the diameter, the stronger the softening. Fig. 0] 
compares phonons for the same (11,11) tube at three dif- 
ferent Tq. It shows that phonon softening increases for 
decreasing To. Note that, as for Figs. 13151 Figs. I6I7I also 
plot only branches corresponding to ZF^^.m) crossing T 
and K. 

Finally, due to the phonon softening induced by the 
KA, we observe that in metallic tubes of any chirality 
the LO branch always crosses the TO, in contrast with 
Ref. Ill, where an anti-crossing of the branches is pre- 
dicted for chiral tubes. 

Fig. [HI plots the LO/TO phonon frequencies for a 
variety of chiral and achiral tubes, both metallic and 
semiconducting, in the diameter range 0.8-2.7 nm at 
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FIG. 8: (color on-line) (a) Frequency of the LO and TO 
modes calculated by EZF-DFT at Tc = 315K. The LO 
phonons are in red; the TO phonons in black. Different sym- 
bols identify metallic and semiconducting (SI and S2) tubes. 
Dashed lines represent the fit of the EZF-DFT data, and cor- 
respond to the plot of Eq. 1281 (b) LO and TO phonons of 
semiconducting tubes only. For S2 tubes TO phonon is ~ 2 
cm~^ higher than in SI tubes. No family dependence is ob- 
served for the LO phonon 
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FIG. 9: (color on-line) Calculated LO and TO modes 
obtained by EZF-DFT for the semiconducting tubes 
(19, 0), (18, 2), (17, 3), (17, 4), (16, 5), (15, 7), (13, 9), and 
(12, 10), listed by increasing chiral angle. All tubes 
have diameter d ~ 1.5 nm. Filled symbols refer to SI 
tubes (mod(n — m, 3) = 1); open symbols to S2 tubes 
(mod(n — m,3) = 2). In S2 tubes the TO mode frequency is 
~ 2cm^^ higher than in SI. 
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FIG. 10: Diameter dependence of LO mode in armchair and 
metallic zigzag SWNTs calculated by EZF-DFT for Tc =4K, 
79K, 315K,3150K. 



E. Analytic results 

So far we relied on numerical DFPT calculations. We 
now show that the key results of the previous Section 
can be derived and explained by using an analytic model. 
From Eqs. 18191 we have: 



(19) 



KAs are due to the presence of non-analytic terms in 
the dynamical matrix and their shape is determined by 
0q, whose expression, within a static approach, is given 
by Eq. [161 ©q can be determined analytically for any q 
at To = OK, or at q — r,K for any To. This provides 
simple formulas for the determination of the metallic- 
semiconducting transition associated with the PD and 
for the KA description in SWNTs. 



1. Kohn Anomalies at T=OK and Peierls Distortions 



To=315K. As expected, we observe no difference in the 
phonon frequencies for TO modes of metallic and semi- 
conducting tubes. Furthermore, confinement does not 
induce any diameter dependence of the frequencies of the 
TO modes. 

Due to the KA presence, LO frequencies are very dif- 
ferent in metallic and semiconducting tubes. In metallic 
tubes there is a strong diameter dependence, with higher 
softening for decreasing diameter. This softening is en- 
tirely due to confinement. Chirality does not affect the 
LO phonons. This is shown in Fig. [9l which plots LO 
and TO modes for all semiconducting tubes with d ~ 1.5 
nm as a function of chirality. 

However, even if there is no dependence on the chi- 
ral angle for the LO modes, a family dependence is ob- 
served for the TO mode in semiconducting tubes. In- 
deed, Fig. [5] (a, b) and Fig. [5] show that TO phonons are 
more scattered than LO. By labeling SI SWNTs with 
mod(7i-m, 3) = 1, and S2 SWNTs with mod(n-r7i, 3) = 
2, it is possible to observe that the TO frequency in SI 
tubes is always ~ 2 cm~^ smaller than in S2 tubes. The 
TO frequency for metallic tubes is in between that of SI 
and S2 tubes. 

Fig. [TO] plots the LO phonon of metallic tubes as a 
function of T^. Since no chirality dependence is observed 
for ro=315 K, only armchair and zigzag tubes are con- 
sidered. The LO softening increases for decreasing Tc, 
and the trend is stronger for smaller diameters. 

We stress again, as in Section II-D, that the temper- 
ature effects calculated here with EZF-DFT are related 
only to the electronic temperature. Therefore, for a di- 
rect comparison with experiments, the data in Fig. [TO] 
need to be corrected with the anharmonic effects. 



To determine the KA shape at zero temperature we 
put To = in Eq. [TO] This gives /k = 1 and /k = for 
states respectively below and above the Fermi energy. 
Eq. [T6]can then be integrated analytically to give: 

c>T,n/K = TTTT. -I- Clo/K' 



^LO/K — 
©TO = 



Ar\/3ag2(i^^)F 



(20) 



where (^'r/x)? defined in Tab. |T] and we assume the 
L and R bands to be linear with slope (3 = 5.25 eV A. 
More details on the integration are given in Appendix B. 
From Eq. [50] we observe that 8to does not depend on q. 
This explains the absence of KA for the TO mode. On 
the other hand, Olo/k has a logarithmic dependence on 
q, which explains the presence of KA in the LO branch. 
Using Eq. [TO] we get: 



"^LO/K 

ar/K 



"r/K 



In 1^1 + Clo/k, 



2V3Ar/Kao(^r/K)F 



7r2/3Af 



(21) 



where q is measured from F for the KA at T and from 
K for the KA at K; ar = 7.89 x 10^ cm-^A"!; q;k = 
7.96 X 10^ cm~^A~^; Clo/k account for all the non- 
divergent terms of the dynamical matrix. 

According to Eq. [HI KA in SWNTs have a logarith- 
mic shape. For the LO and K modes at Te = OK, lo^ 
is negative, which gives imaginary phonon frequencies. 
This means that the lattice undergoes a permanent dis- 
tortion, i.e. a Peierls distortion. 

Fig. [H plots Eq. [Ufor the KA at = OK and com- 
pares it with EZF-DFT calculations in the limit Tc ^ OK. 
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FIG. 11: KA at T and K for a (11,11) tube. Triangles cor- 
respond to EZF-DFT at rc=4K. The solid line is the plot of 
Eq.[2T] 



Using Eq. (22 and fitting Clo/k to the EZF-DFT data, 
we get: 

ujIo = 5.3 X lO'^lnq + 2.70 X 10^cnl"^ 

= 5.3 X 10'^ln(7 + 1.94 X lO^cm^^ (22) 

From Fig. [Til Eq. perfectly reproduces the phonon 
dispersions obtained by EZF-DFT. 



2. Electronic Temperature dependence of Kohn anomalies 
and Peierls distortion temperature 



We now fix q = r or q = K and study the dependence 
of the phonon frequencies on Te, which changes the oc- 
cupation functions /q,n- For q — r,K, assuming linear 
L and R bands with slope (3 = ±5.52 eVA, remembering 
that K -I- K = K', and setting ep = 0, we have: 

e(K-i-k'),L = £(K-i-k'+r),L = e(K+k'+K),L = k'/3 (23) 
e(K-i-k'),-R e(K+k'+r),_R, = e(K+k'+K),fl = -k'^ (24) 

The occupation factors /q_„ depend on q and n only 
through e„(q) — ep. Defining x = I3k' /Tk-Q, and 9{x) = 



To (K) CC (eV^) CC (cm"^ 



9612 
2646 



0.039 
0.027 



2.54 X 10" 
1.76 X lO'' 



TABLE IV: values of To and CCr,K obtained fitting the EZF- 
DFT data for tubes calculated at T^ = 4K, Tc = 77K Tc = 
315K Te = 3150K and for d in the range 0.8-2.4 nm with 
Eq.m 



1/(1 -I- e^), Eq. [TBI becomes: 

^r/K2V3a§(I?2^K>F 



e 



r/K 



ah 

e(x) 



-dx, 
(25) 



■ Tk-p 



where the bands' symmetry allows to replace J ^''Jj with 
2 //'^ . This gives: 



"r.K , T 
— In — 

d To 



r,K, 



(26) 



where ctq is the same as in Eql2H To is the electronic 
temperature for which the contributions of Oq are null, 
C'Cr.K is the value of Wp^j, at Te = To, and contains 
all the contributions from Q^. Details on how Eq. [^D 
is derived from Eq. [16| are given in Appendix C. Eq. [26] 
can be used to fit simultaneously all the phonon frequen- 
cies calculated by EZF-DFT as a function of To. The fit 
parameters are listed in Tab. IIVI 

The electronic temperature, TpD, for which a SWNT 
undergoes a Peierls distortion, can be obtained by setting 
cjq=0 in Eq. [211 which leads to: 



TpD = Toe °r,K 



(27) 



As an example. Eg . [27l apphed to the (6, 6) tube {d = 
0.8nm) gives Tpo = 6.08 x lO'^K for the phonon at T, 
and TpD = 0.5 x 10~^K for the phonon at K. Tpd expo- 
nentially decreases with the tube diameter. This implies 
that, unlike the ultra-small tubes studied in Refs. [T9ll20l . 
for the SWNTs generally used in experiments, TpD is 
always smaller than 10^^ K. 

Eq. [211 defines the diameter dependence of the LO 
phonon frequency of metallic tubes (^lq) ^t any To, and 
can be used to fit the EZF-DFT data of Fig. M We 
observe that all the EZF-DFT data in Fig. [8] can be rep- 
resented by very simple functions of the tube diameter. 
Indeed, neglecting the very small difference observed for 
the SI and the S2 tubes, the frequency of the TO mode of 
both metallic and semiconducting tubes (wto) is diam- 
eter independent, while the calculated frequencies of the 
LO mode for semiconducting tubes (wlq) ^re inversely 
proportional to the tube diameter. This can be summa- 
rized by the following equations: 



d ■ 



Vs 



d ' 



WTO = t (28) 
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where z/m = —77.33 cm"-'^ nm, rju = 1597 cm~^, i^s = 
13.78 cm~^ nm, and rjs = ^ = 1579 cm~^, which per- 
fectly reproduces the phonon frequency we calculate for 
the r — E2g mode of graphene. The functional form of 
uj^Q has been obtained from Eq. [26] using the relation 
\/l + X ^ 1 + x/2. The expression for is thus com- 
pletely analytical, but is valid only for Te = 315i4r. On 
the other hand, the expressions for ujto and w^o ^-re em- 
pirical, but are valid for any T^. Eqs [28] are plotted in 

Fig. m 



V. CONFINEMENT EFFECTS WITHIN 
DYNAMIC DFPT 

In Section II-E, we argued the possible presence of dy- 
namic effects for one-dimensional metallic systems. To 
the best of our knowledge, the dynamic nature of phonons 
has always been neglected in the SWNT literature, until 
we recently pointed this out^^. In this Section, we show 
that dynamic effects are present in metallic SWNTs, and 
that they induce significant changes to the KA occur- 
rence and shape. Therefore, they must be included to 
explain the experimental data. In fact, SWNTs are one 
of the first real materials for which is a significant differ- 
ence between static and dynamic DFPT is detected. 

If the static approximation is relaxed, the dynamical 
matrix is described by Eq. [H and its non analytic part 
by Eq. [TO] Assuming T = 0, Eq. [TO] can be solved by 
using a procedure similar to that used to obtain Eq. 1201 
For the LO phonon, assuming fc ^ |g ± ^^|, we get: 



0.00 0.01 0.02 0.03 0.04 0.05 0.06 0.07 



while for the TO: 



©TO — — 



Arpq^Vial2{Dl,)F 



(29) 



(30) 



where Wlo/to is the phonon frequency of a phonon of 
wave vector g ~ F belonging to the same branch of the 
LO/TO phonon at q = 0. These equations mark a major 
departure from the previous results. Qq in the dynamic 
approach is very different from the static case. Note that 
the static case (Eq. I^Hl) is immediately recovered setting 
ftwLO/TO = in Eas. l29l30l 

The dynamic effects are qualitatively very different for 
LO and TO phonons. Eq. [29] diverges for q = hcUq/P, 
and not for g = as in the static case. This gives the 
first major result that the KA for the LO mode is not at 
F, contrary to what is predicted by any static approach, 
but a,t q = ±huJLo/P- 

The dynamic effects are even bigger for the TO 
phonons. Eq. [30] shows that 8to is zero for q ~ 0, 
but diverges to -l-oo for q [hw-io/P)^ and to —oo for 
q — > [hiOTo / (3)'^ ■ This is quite different from the static 
case, where Oto is a negative constant. This results in 
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FIG. 12: (color on-line) KA for the LO and TO branches of an 
(11,11) tube close to T. Phonon frequencies are calculated by 
means of static EZF-DFT (open dots), or from the numerical 
integration of the dynamical matrix in its static (dashed line) 
and dynamic (solid line) expression. To = 315K. 



a significant TO upshift with respect to the static case. 
This can be described by a simple equation, the deriva- 
tion of which is given in Appendix D: 



dyn 



25 



(31) 



where the frequency difference is in cm~^ and d in nm. 

To summarize, the static approach; (i) fails to describe 
the position of the KA for the LO mode, (ii) underes- 
timates the TO phonon frequency, (iii) misses the TO 
Kohn anomaly at q = Hloto I P- 

The phonon dispersion of metallic SWNTs close to the 
KAs is obtained by adding the contributions from the 
analytic and the non-analytic part of the dynamical ma- 
trix. The contribution of the non-analytic part is given 
by the numerical integration of Eq. [TO] Due to the pres- 
ence of the phonon energy in the expression of Qq, Ea.[TOl 
has to be solved self-consistently. The contribution from 
is easily obtained from EZF-DFT. The definitions of 
depend on the value chosen for k. In our 
calculations, we chose k such that = 1.0 eV. 

Fig. [TO] compares the phonon frequencies of the LO and 
TO branches of an (11,11) tube calculated using EZF- 
DFT with those obtained from the numerical integration 
of Eqs. I16I10I Fig. [TO] shows major differences between 
the dynamic and the static approach in the region of the 
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FIG. 13: (Color on-line) KA of the LO and TO branches for 
tubes of different diameters calculated at Te = 315K including 
dynamical effects. 



0.00 0.01 0.02 0.03 




0.00 0.01 0.02 0.03 



Phonon Wavevector (27i:/ajj units) 

FIG. 14: (color on-line) KA of LO and TO branches of tubes 
with d — 1.5 nm calculated for To ranging from 30 lOOOK, 
including dynamic effects. 



A. Dynamic effects in graphene and graphite 



phonon dispersions affected by KA. As expected, in the 
dynamic case, the KA in the LO branch is shifted from 
g = to a finite wave- vector. This causes a significant 
change close to F. However, Fig. [12] also indicates that 
these changes are confined to a very small BZ region. In- 
deed, for q > 0.02 in 27r/ao units, the static and dynamic 
results are indistinguishable. 

It is important to observe that Eq. [271 which was de- 
rived within the static approximation, is also valid in the 
dynamic case. Indeed, for Tc 0, metallic SWNTs un- 
dergo a PD, meaning that the frequency of the phonon 
affected by the KA goes to zero. For fioiLo/TO = the 
dynamic equations reduces to the static case, so TpD is 
the same. 

Figs. I13I14I plot respectively the KA shape in the LO 
and TO branches as a function of diameter and To. The 
general trends observed within the static approximation 
are preserved also after the introduction of the dynamic 
effects. In particular, the LO phonon softening still in- 
creases for decreasing diameter and Tq. In addition, the 
phonon softening/hardening of the TO branch increases 
for decreasing diameters, but has no dependence on Tc. 



In this Section we have shown that dynamic effects 
induce huge modifications in the phonon dispersion of 
metallic nanotubes close to the Kohn anomaly. The 
demonstration of the importance of dynamic effects in 
one-dimensional systems opens the discussion on whether 
similar effects can be relevant in materials of higher di- 
mensionality, graphite and graphene in particular. In- 
deed, it is possible to prove that dynamic effects can be 
crucial also for the correct description of the vibrational 
properties of bi-dimensional metallic systems. We have 
recently demonstrated, both theoretically and experi- 
mentally, that for electrically doped graphene the Born- 
Oppenheimer approximation spectacularly fails^?'^^. In 
particular the evolution of the G peak position as a func- 
tion of gate-bias can be explained only by taking into 
account non-adiabatic effects, as will be reported in de- 
tails elsewher e 



VI. CURVATURE EFFECTS 

Curvature effects encompass all the differences between 
graphene and SWNTs due to the geometrical distortion 
of the C-C bonds in the tubes (Eq. [T^ . In common with 
quantum confinement, curvature also splits the double 
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degenerate T — E2g phonon of graphene into two distinct 
phonons in SWNTs. This split occurs because the differ- 
ent strengths of the chemical bonds along the tube axis 
and the tube circumference result in different phonon 
frequencies for modes polarized along these two direc- 
tions^!^. The folding of graphene to give the SWNT 
cylindrical shape results in: (i) change in atoms sepa- 
ration (ii) loss of C-C bond planarity, (iii) mixing of a 
and TT states, giving a sp^ / sp^ character to the chemical 
bonds. 

When the purely sp^ bonds of graphene are deformed, 
they assume a hybrid sp'^ /sp^ character. Since sp^ bonds 
are stifFer than sp^ bonds, mixing the orbitals results 
in bond softening. This mixing is proportional to cur- 
vature, so it is minimum for C-C bonds parallel to the 
axis, and maximum for bonds oriented along the circum- 
ference^i^. Since longitudinal modes deform C-C bonds 
parallel to the axis, while transverse modes deform C-C 
bonds along the circumference, one expects curvature to 
soften the modes polarized along the tube circumference 
more than those polarized along the axis. This softening 
should increase with the sp'^ /sp^ mixing, and, thus with 
the SWNT curvature. 

In this Section, we investigate the effects of pure cur- 
vature on the LO-TO splitting. In Eq. [12] the dynamical 
matrix of a SWNT (G)t) is written as the sum of the 
dynamical matrix of graphene (6g) and two additional 
terms describing the effects of curvature (Ocurv) and con- 
finement (Oconf). EZF-DFT describes the contributions 
deriving from graphene and from the confinement effects 

(efl-t = eG + eco„f), thus: 

Ocurv = Gt - Og - Oconf = Gt - Ot"'- (32) 

According to Eq. [321 the effects of pure curvature on a 
given SWNT can be obtained as the difference between 
calculations performed on the actual SWNT and the re- 
sults obtained by means of EZF-DFT. 

The frequency dependence on the diameter for the 
modes derived from the graphene T — E2g was calcu- 
lated by several authorsiiii^i^l. As already pointed out 
in Section IV, previous calculations of the LO phonons 
in metallic SWNTs are either not converged with respect 
to the k-point sampling or use a fictitious electronic 
temperatureiiiiSiiS instead of the real one, and thus do 
not properly describe the KA effects. However, for semi- 
conducting SWNTSjKA are not present. Thus, the DFT 
results of Refs. llTlllSl for LO and TO modes in semicon- 
ducting SWNTs are reliable, and can be used to estimate 
the curvature effects for LO and TO modes. 

To verify the equivalence of the calculations of Refs. 
[TtIITsI with the present, the phonon dispersions of 
graphene of Ref. [3 are compared with those of Ref. [l3 
(which is obtained with the same method and code as 
the present calculations). Away from the KA, the only 
noticeable difference is a 16 cm~^ rigid upshift of the 
highest optical branches. Thus, the two calculations are 
equivalent (after a downshift of 16 cm~^ ) and is possible 



to extract the effects of pure curvature by comparing the 
results of Refs. [l3ll with the present EZF-DFT calcu- 
lations. 

Assuming the curvature effects as a perturbation, we 
define 

Aw = tJT - Wfiat, (33) 

where wt is the phonon frequency calculated on an ac- 
tual tube and Wflat is the EZF-DFT result. In FigfTSl we 
plot Aw as the difference between the phonon frequencies 
of ReffTs' and the Wflat from EZF-DFT as represented 
by Eq. [28] We consider only those modes which do 
not present any dependence on To. In particular, for the 
LO mode, we consider only semiconducting tubes, while 
for the TO mode, we consider both metallic and semi- 
conducting tubes. For these modes, it is not necessary to 
distinguish between static and dynamic EZF-DFT. From 
a fit to the points of Fig. [T51 

Aa;To(rf) = 

AwLo(d) - (34) 

where Qro = 25.16 cm^^ nm^, Clo — 12.0 cm^^ nm^, 
and d is the tube diameter in nm. We assume that Eq. (341 
applies also to the LO mode of metallic tubes. 

Fig. [TBI plots the diameter dependence of the LO/TO 
phonons as predicted considering only the effects of con- 
finement (static and dynamic) or after the introduction 
of the curvature effects (added using Eq. [34]). In metalhc 
tubes, the relaxation of the static approximation intro- 
duces a diameter dependence for the TO phonon, and 
reduces the softening of the LO mode. This is a further 
confirmation that in SWNTs dynamic effects cannot be 
neglected. The most notable modification introduced by 
curvature is the strong diameter dependence of the TO 
mode in semiconducting tubes, which, according to pure 
confinement, was predicted to be diameter independent. 
Interesting, curvature almost perfectly compensates the 
effects of confinement for the TO mode in metallic tubes 
and for the LO mode in semiconducting tubes, resulting 
in almost diameter-independent phonons. Finally, the ef- 
fects of curvature in the LO modes of metallic tubes do 
not substantially modify the trend due to pure confine- 
ment, and result only in a correction of the phonon fre- 
quencies. Fig. [T6]is extremely important, and shows that 
curvature has the effect of inducing a noticeable phonon 
softening in small diameter SWNTs, for both LO and TO 
modes. 



A. LO-TO splitting 

In Sections IV and V, we observed that quantum con- 
finement has the effect of splitting the graphene T — E2g 
phonon into a LO and a TO mode. Literature reports 
indicate that in chiral tubes curvature effects may induce 
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FIG. 15: (color on-line) Curvature effects on LO and TO 
phonons. Dots re pre sent the difference between the phonon 
frequencies of RefEi and EZF-DFT of Section IV. (Ea l28)l . 
The lines are fit to the data done using Eg 1341 
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FIG. 16: (color on-line) Comparison between LO and TO 
phonons. (a) for metallic tubes the phonon frequencies 
are calculated using EZF-DFT within the static approxima- 
tion (dotted line), EZF-DFT including the dynamic effects 
(dashed line), and dynamic EZF-DFT correct with the cur- 
vature effects (solid lines), (b) for semiconducting tubes it is 
not necessary to distinguish between static and dynamic cal- 
culations; dotted lines represent the EZF-DFT calculations, 
solid lines represent EZF-DFT data corrected for the effects 
of curvature. 



FIG. 17: (color on-line) Pattern of atomic displacements as- 
sociated to a phonon derived from the graphene T — E2g mode 
(a) as predicted by EZF-DFT and (b) as calculated in Ref . [Sll. 



deviations from this behavior-i^i^. In Ref. |5lj, it was ob- 
served that the eigenvectors of the LO and TO phonons 
in the (9,3) tube are aligned to the C-C bonds, and devi- 
ate by 16° from the axial and circumferential directions. 
However, in the same paper, it was also reported that 
for the (8,4) tube the deviation is reduced to 2°. Cal- 
culations done in Ref. Is' show that for the (12, 6) tube 
the deviation from purely LO and TO eigenvectors is 4°. 
Using a symmetry adapted TB scheme, calculations per- 
formed in Ref. [3l| on 300 different SWNTs indicated no 
substantial deviation of the atomic displacements from 
the axial and circumferential directions. Fig. [T7] shows 
the very different predictions given by our EZF-DFT and 
the calculations of Ref. [5l|. 

The orientation of the atomic displacements in 
phonons derived from the T ~ E2g of graphene can be 
understood by considering the combined effects of con- 
finement and curvature. The frequencies of the T — E2g 
modes of graphene are degenerate, and the correspond- 
ing eigenvectors can be oriented along any couple of in- 
plane directions. Curvature and confinement both have 
the effect of removing the degeneracy. The effect that 
curvature and confinement have in orienting the atomic 
displacements is proportional to the effect they have in 
determining the LO-TO splitting. These contributions 
can be read from Fig. 1161 by comparing the curves be- 
fore and after the introduction of the curvature correc- 
tions. In metallic tubes the effects of confinement are 
dominant over curvature, thus the modes derived from 
graphene T — E2g are expected to be almost perfectly 
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LO and TO. In Section IV, and in Tab. IIIII we pointed 
out that the calculations in Ref. [5l| were done using a 
very limited number of k-points, and thus failed to de- 
scribe the effects of confinement. This is also supported 
by the results presented in Ref. 18. Indeed, calculations 
of Ref. H are converged with respect to the number of 
k-points, and show that for the (12,6) tube the modes 
are almost perfectly LO and TO. 

For semiconducting tubes, the effects of curvature are 
dominant over the effects of confinement. In this case, 
the almost perfect LO and TO character imposed by 
confinement can be substantially perturbed by curva- 
ture. However, it is important to remember that the 
splitting caused by curvature is due to the different force 
constants along the tubes axis and circumference. As a 
consequence, phonons that present a strong curvature in- 
duced splitting have to be oriented along directions for 
which the difference between the force constants are sig- 
nificantly strong. This means that the atomic displace- 
ments have to be substantially aligned along the tube axis 
and the tube circumference. This suggests that the devia- 
tions from purely axial and circumferential displacement 
have to be modest, and that it should be always possible 
to classify the modes as almost LO and almost TO. 



VII. INTERPRETATION OF THE RAMAN G 
BAND IN NANOTUBES 

We now show that combining curvature, confinement 
and dynamic effects we can interpret the SWNT G peak. 

The Raman spectrum of SWNT shows strong fea- 
tures in the 1540-1600 cm~^ spectral region^i^i^iSOiSi^^ 
This is the same region of the G peak in in 
graphenei^, graphite^ and amorphous and nanostruc- 
tured carbons^^iS^. 

As shown in Fig.[lHl the G band of SWNTs consists of 
two main peaks: G"*" and G^ In semiconducting tubes 
both peaks are sharp, while in metallic tubes the G^ 
is broader and downshifted^'^'^^'28^'20^^^^^^. 
For both metallic and semiconducting tubes, the position 
of the G~ peak depends on the tube diameter, having 
lower frequency for smaller diameters, while the posi- 
tion of the G"*" peak is substantially diameter indepen- 
dentS^. As discussed in Section VI, curvature affects 
more strongly the circumferential modes, thus the G^ 
and G~ are commonly assigned to LO (axial) and TO 
(circumferential) modes, respectively i^'^^i^^'^^i^^'^^ 

Conflicting reports exist on the presence and rela- 
tive intensity of the G~ band in isolated versus bundled 
metallic tubes. Some groups report that the intensity 
of this peak in isolated tubes is the same as in bun- 
dleSj^iS^iS^ while others observe that it is smalleriSLS^iSi 
or even absent^. The downshift and the broadening of 
the G~ peak in metallic tubes is commonly attributed 
to the onset of a Fano resonance between plasmons and 
the TO phononi22ii21i22ii2ii2^. Such phonon-plasmon cou- 
pling would either needSS, or not need?2,'M a finite phonon 



wavevector for its activation. The theory of Refs. [qoIIq^ 
predicts the phonon-plasmon peak to be intrinsic in sin- 
gle SWNT, in contrast with Ref. [H On the other 
hand, the theory in Ref. fol requires several tubes (>20) 
in a bundle in order to observe a significant G~ inten- 
sity, in contrast with the experimental observation that 
bundles with very few metallic tubes show a significant 
Ref. M also predicts a G" upshift with 
number of tubes in the bundle, in contrast with Ref. 113, 
which shows a downshift, and with Refs. [sgIIooIIqsI . which 
show that the G~ position depends on the tube diameter 
and not bundle size. Finally, the G~ position predicted 

by Refs.[9l9i is at least 200 cm lower than that mea- 
sured 5.6 .86.87.88.89.90.91.93.94.95.96 

Thus, all the proposed theories for phonon-plasmon 
coupling22i22ii2S are qualitative, require the guess of sev- 
eral quantities, and fail to predict in a precise, quanti- 
tative, parameter-free way the observed line-shapes and 
their diameter dependence. 

Phonon symmetry proves that there are only 6 Raman- 
active mode at ri^iSO^^ioo. 2 A, 2 Ei, and 2 E2. The 
A modes are totally symmetric, have the highest EPC, 
and correspond to the LO and TO modes we examined 
in the previous Sections. In Figliniwe compare the ex- 
perimentally measured G+ and G^ peaks to the calcu- 
lated frequencies of the LO and TO phonons. The plot 
of the phonon frequencies corresponds to the results of 
static and dynamic EZF-DFT plus the curvature correc- 
tion Auj. 

For semiconducting tubes, our calculations support the 
G+, G^ assignment to the LO and TO phonons. 

For metallic tubes, our LO and TO positions reproduce 
the diameter dependence of the G~ and the G+ peaks 
respectively. Our assignment of the G peaks in metal- 
lic tubes is thus the opposite of that in semiconducting. 
This is emphasized in Fig |181 where we represent the cor- 
respondence between the Raman peaks and the LO/TO 
phonons. 

Our calculations prove that the G~ softening can be 
precisely predicted by simply considering the coupling of 
the LO phonon with single-particle electronic excitations. 

Furthermore, in Ref|70|, we have shown that our DFT 
calculated EPC also account for the measured FWHM. 
Vice- versa the experimental FWHM can be used to mea- 
sure the EPCffl. It is important to remind here that 
only single resonance Raman can detect a G~ with large 
FWHM in metallic nanotubes---. The same metallic nan- 
otubes measured in double resonance would have a sharp 

g - 70.93 

It is also interesting to consider if Raman scattering 
can allow to probe regions very close to the anomalies in 
FigmmSmH This is not possible in the single resonance 
approach, since it only probes q = 0. This might be pos- 
sible in double resonant Raman scattering. In this case 
the q = rule is relaxed, and the scattering process in- 
volves phonons with q 7^ q 49,ioi ^ However, only phonons 
with q > hujr / P satisfy the conservation of energy and 
momentum in a double resonant process. Fig. [T^] shows 
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FIG. 18: Raman and G~ peaks for semiconducting and 
Metallic SWNTs (spectra adapted from Ref. HI). In semi- 
conducting tubes the two peaks have a Lorentzian shape and 
a FWHM of ~ 12 cm~^. In metallic tubes the G~ peak is 
downshifted and much broader. Our assignment of the spec- 
tral features is also indicated. 



that for such q the anomaly is missed and the predicted 
peak position is very similar to the q — case. However, 
crossing from single to double resonant would result in 
a significant downshift of the TO mode, as for Eq|3T] 
This is consistent with the experimental data presented 
in Ref. ^ 

Further validation of our calculations with experiments 
is provided by the analysis of the temperature depen- 
dence of the G"*" and G~ positions. 

As stressed in the previous Sections, here we do not de- 
scribe the frequency shift induced by the presence of an- 
harmonic effects. However, the effects of anharmonicity 
are expected to be the same for the LO and TO phonons. 
Thus anharmonicity does not affect the relative position 
of these two phonons, which can then be altered only 
by effects related to T^. The temperature dependence 
of the G+ — splitting can be investigated using our 
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FIG. 19: comparison between LO and TO modes of (a) metal- 
hc and (b) semiconducting tubes and the experimental posi- 
tion of the G"*" and G~ peaks of the Raman spectra of SWNTs. 
The position of the LO and TO modes is calculated using 
EZF-DFT and adding the curvature correction. In panel (a) 
continuous lines represent the results of dynamic EZF-DFT 
calculations, while dotted lines refer to static EZF. Open and 
filled symbols represent respectively the measured position of 
the G"*" and G~ peaks. Data from Oron et al.-^^ (squares), 
Jorio et al.** (circles), Maultzsch et al.i^, (triangles), Doorn 
et al.-SS (diamonds), Meyier et al.— (hexagons), and Brown at 
al.— (pentagons) 



analytical model. 

The difference between the frequency of the LO and 
TO modes for metallic and semiconducting tubes of sim- 
ilar diameter d = 1 nm is plotted as a function of Tq in 
FigEOl For the semiconducting tube no effects related 
to To are expected, thus the splitting is temperature- 
independent. For metallic tubes, the static and the dy- 
namic model give very different predictions. According to 
the static model, the frequency split should have a strong 
logarithmic dependence on temperature, Eq. 1261 On the 
contrary the dynamic model, obtained by the numerical 
integration of Eq. [TOl predicts that the splitting (i) has a 
maximum at ~ 500K, (ii) it saturates for — > 0, (iii) 
for To 00 it decrease monotonically, and (iv) in the 
range T = ^ lOOOK it varies less than 10 cm'^ This 
variation is function of the tube diameter, and is plotted 
in Fig. [21] We observe that this value is extremely small 
(~ 2 cm^^) for large diameters, but rapidly increases for 
d < 1.5 nm. 
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FIG. 20: Tc dependence of LO-TO splitting for a metallic 
and a semiconducting tube, calculated using EZF-DFT and 
including curvature corrections. The solid line corresponds 
to results of dynamic calculations on a metallic tube; the 
dashed line correspond to static calculations. The splitting 
for a semiconducting tube is represented by the dotted line. 



Experimental data from Refl77l. obtained from tubes 
having d ^ 1.0 nm, are in excellent agreement virith the 
prediction of the dynamical model. This confirms the 
validity of our model and the inadequacy of a static ap- 
proach for the investigation of the vibrational properties 
of SWNTs. 



VIII. CONCLUSIONS 

We presented a detailed theoretical investigation of the 
optical phonons of SWNTs. We first analyzed the static 
DFPT approach and singled out the effects of curva- 
ture and confinement. Confinement plays a major role 
in shaping SWNT phonons and is often more relevant 
than curvature. We presented an electronic zone fold- 
ing method allowing the evaluation of confinement effects 
on phonon-dispersions and electron-phonon coupling of 
SWNTs of any diameter and chirality, and for any elec- 
tronic temperature. We investigated Kohn anomalies and 
Peierls distortions in metallic SWNTs with diameters up 
to 3 nm and in a Te range from 4 to 3000 K. We then pre- 
sented a simple analytic model exactly accounting for all 
the static DFPT results. Finally, we proved that non adi- 
abatic dynamic effects, beyond the Born-Oppenheimer 
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FIG. 21: difference between the saturation value at Tc = 
OK and the maximum value at Te ~ 500K of the frequency 
difference of the G^ and the G~ Raman peaks as a function of 
the tubes diameter. The values are calculated using dynamic 
EZF-DFT and including curvature corrections. 



approximation, induce significant changes in the occur- 
rence and shape of Kohn anomalies. We have shown 
that it is necessary to consider dynamic effects in order 
to correctly describe the phonon dispersions of SWNTs, 
in contrast to what happens in most other materials. Fi- 
nally we discussed the interpretation of the Raman G 
peak in SWNTs. Only by combining dynamic, curvature 
and confinement effects we can reach an agreement be- 
tween experiments and theory. In metallic SWNTs, we 
assign the G"*" and G~ peaks to the TO (circumferential) 
and LO (axial) mode, the opposite of semiconducting 
SWNTs. 
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X. APPENDIX A 

Here we derive the the relation between the SWNT 
and graphene EPC presented in Section III, Eq. [T31 
By definition of the EPC (Eq. O: 

dH 

D{k+q)nMm — (l/'k+q.n | ^^Eq.r) I V'k.m) , (35) 

where H is the system Hamiltonian, R represents the 
atomic coordinates, {ip^.n) is the electronic state of 
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wavevector k and band n, and Sq is the polarization 
of a phonon of wavevector q and branch 77. We define 
\ipi^^n) 3-nd £q ,j the wavefunction and polarization nor- 
malized to 1 on the unit cell. We also define IV'k'n) ^^'^ 
the wavefunctions and the polarizations normalized 
to 1 on a supercell composed by N unit cells: 

1 



V i V 



-q,r;- 



The EPC calculated on the supercell is then: 



^(k+q)ri,k,m — \rk+q,nl ^jj^ ^q,J) I ^k,7i 



(36) 



(37) 



This integral is equivalent to the sum of N unit cell in- 
tegrals: 



^(k+q)n,k,m 



N{ 



k.rn V 



^(k+q)ri,k,m- 



(38) 



Considering a SWNT with a unit-cell containing 
graphene unit cells, Eq. [T3]is immediately obtained. 

Eq. [37] states that the EPC depends on the supercell 
choice, as a consequence the EPC per se is not a physical 
observable. In general, a physical observable is given by 
the product of the square of the EPC times an electronic 
density of states p: 



P(e) 



1 



(39) 



where the sum is performed on N]^ wavevectors k and 
on the bands n. e„(k) is the energy of \ilj\i,n)- One can 
easily check that p is proportional to TV, the size of the 
super-cell. Thus, the product \D\'^p does not depend on 
the super-cell chosen and can correspond to a physical 
observable. 



XI. APPENDIX B 

We evaluate Eq. [11] to obtain the shape of the KA 
at T = 0, presented in Section IV-E. For T = OK, the 
occupation factors / in Eq. [TH] are or 1 . For q F the 

non-analytic component of the dynamical matrix is then: 



e 



q~r 



At Vial 

TT^d 



^ l^(K+k'+q)7r*,(K+k')xP 
k e(K+k')7r - e(K+k'+q)7r' 



dk', (40) 



where tt and tt* indicate bands occupied or empty, respec- 
tively (Fig. [2]). Neglecting the dependence of \D\ on k' 
and assuming a linear dispersion for the electronic bands 



e(K+k')7r — 

0q~r 



1 



Pk' and e(K+k')7r* = +(3k', Eq.[40l becomes: 
- \Dll\' 



-k P{2k' + q) 



dk' 



-9 

+k 



-(3q 



-dk' 



\D 



RL\ 



1 



-(3{2k' + q) 



dk' 



(41) 



Here, the subscripts of D indicate whether we are consid- 
ering electronic bands L or R (Fig. The integrals can 
be easily evaluated. Considering In |2fc+9| ^ In |2fc| = Cj. 
and \Dlr\^ = \Drl\'^, we have: 



e 



q~r 



n^df] 



\DLR\'Mq\^C-,)-\DLL\ 



(42) 



We now distinguish between the two LO and TO 
branches at q P. From Tab.|Il we get |-Dl_r| = 2(£>p)f 
for the LO branch and for the TO, while I-DllI = 
2{DI)f for the TO and for the LO, thus 



Olo 
9to 



ArViap{Dl)F 



ArV3ap{Dl)F 

7T^d/3 



(43) 



where Olo ( Oto ) is the non analytic component of 
the dynamical matrix related to the LO (TO) branch. 
C'^ — Cj,ArV3al2(D^)p/'K'^dp contains all k dependent 

terms. For g — > Olo ~^ —00, while ©to is a constant 
independent on q. This explains the occurrence of the 
KA only for the LO branch and its absence for the TO 
branch in static DEPT. 

The phonon frequencies uj are obtained by = 
3?e{(9 4- Q'^'^)/M}, where M is the carbon mass and the 
analytic component of the dynamical matrix Q^"^ is ex- 
pected to have a very weak dependence on q. Conclud- 
ing, at T=0, for q '--^ F, the frequencies of the LO and 
TO branches are 



ArVSal2{Dl)F 



■'LO 



2 

WrpO 



-In I, 



TT^dPM 

Clo, 
ArV3al2{Dl)F 

TT^dPM 



In \q\ + Clo 



C' 



TO 



(44) 



(45) 



where the constants Clo /to include all the contributions 
from the analytic part of the dynamical matrix O''" and 
all the non-divergent terms of 0. For the numerical eval- 
uation of ar we use: Ajp — 2, qq ~ 2.46 A, /? = 5.52 
eV A, {DI)f = 45.60 eV2A-2, M = 12.01 a.m.u. The 
corresponding equations for the KA at K are obtained in 
a completely analogous way. 
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XII. APPENDIX C 

We determine the dependence on of the KA phonon 
frequencies at q = F and q = K. We define: 



where In(xo) = - F{Q) + In(ii) and Tq = 

Considering that the analytic part of the dynamical ma- 
trix (Eq. [19]) as independent from T, Eq. [26] is then ob- 
tained. 



/(^) = ^ 
Using Eq. [ST] Eq. [25] becomes 

<Jr/K - — 1— 



f{x)dx, 



(46) 



(47) 



where x = f3k/{Tk^). We observe that lima;^o/(a^) = 
and that for x — > oo /(x) ~ — Since fc is chosen in 
the linear range of the tt and tt* bands, /3fc' is a frac- 
tion of eV. Moreover, fee = 8.617 x IQ-^eV/K, so the 
condition j3k' ^ fceT, corresponding to the x — > oo, is 
reached for T up to thousands K, well above any realistic 
experimental condition. In Eq. [47]we replace /(x) with: 



XIII. APPENDIX D 



We evaluate the frequency upshift of the TO mode 
due to dynamic effects. From Eq. [3D] we observe that, 
for q=0, -Dto = 0. Thus, from Eq. [19] one obtains that 
for q=0: 



M ' 



TO I 



(51) 



Eg. [171 becomes: 

"r/K 



/(x) for X < xi 
for X > xi 



6 



r/K 



/(x)dx 



Defining F{x) = J f{x)dx, we get: 
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«r/K 
d 

ar/K 



c? X 
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(48) 



(49) 



Xi 



(50) 



Since the contribution of the phonon energy in the de- 
nominator of 0q is negligible for k away from the Fermi 
surface, I?™ is the same in the static and in the dynamic 
case. Using Eqs. [T9l [20] and [211 we get: 



Eqs. [STlandl 



^TO 
1 • ^ 

J gives: 



/, ,stat\2 
I'^TO ) 



ttr 
~d' 
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(52) 



(53) 



Eq. [31]is finally obtained by using \/l~-f 
substituting the numerical value of ar . 



1 -I- x/2 and 
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